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We analyze the perturbative series of the Keldysh-type sigma-model proposed recently for de- 
scribing the quantum mechanics with time-dependent Hamiltonians from the unitary Wigner-Dyson 
random-matrix ensemble. We observe that vertices of orders higher than four cancel, which allows 
us to reduce the calculation of the energy-diffusion constant to that in a special kind of the matrix (f>* 
model. We further verify that the perturbative four-loop correction to the energy-diffusion constant 
in the high-velocity limit cancels, in agreement with the conjecture of one of the authors. 

I. INTRODUCTION 

While spectral properties of static random-matrix Hamiltonians have been studied in detail by various methodsi*^, 
much less is known about quantum-mechanical evolution with a time- dependent random-matrix Hamiltonian: 

i-^if^nm. (1) 

Such a system with the Hamiltonian H{t) belonging to one of the three Wigner-Dyson random-matrix ensembles 
(unitary, orthogonal or symplectic) has been studied by Wilkinson in Ref. 0. The trajectory of the Hamiltonian H{t) 
in the space of Hermitian matrices is assumed to be nearly linear on the time scales relevant for the problem. This 
assumption is justified in the limit of large matrix dimension N: in this limit the energy level spacing A is small, and 
a small variation of the Hamiltonian matrix elements (of order A, which is much smaller than the matrix elements 
themselves) already shifts the energy levels by the order of the level spacing and thus changes the level correlations 
completely. Therefore the relevant time scales are small in iV, and any trajectory H{t) with smooth (independent 
of N) time dependence may be approximated as linear. This is the usual reasoning in the studies of parametric 
level statistics^ which deduces that under such an assumption the spectral correlations acquire universal properties 
(independent of the particular choice of the trajectory). 

Following the traditional notation and for the future possibility of describing different time dependencies of H{t), 
we introduce the time dependence of the Hamiltonian in two steps. First we define the class of linear trajectories 
H{ip) and then let the parameter iphe a given function of time t. In the present paper we restrict our discussion to the 
Gaussian Unitary Ensemble, for which the linear trajectories H(ip) may be defined by the pair-correlation functions 
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Hij{ip)Hki{(p) = — —5a5jk, (2a) 

{H{^)-H{^'))^^{H{v)-H{v'))^^ = 5aSM^ ~ ip'fA^C{Q) + O (^^^^) , (2b) 

where A is the mean level spacing in the center of the Wigner semicircle, and C(0) is the conventional notation for 
the sensitivity of the energy spectrum on the parameter (see, e.g., Ref. 151). 

Two possibilities of the motion along the trajectory H{ip) are of principal importance: 

(i) Linear time dependence ip{t) ~ vt. This is the problem studied by Wilkinson** and also the situation considered 
in the present paper (except in the Section IVHIl where more general time dependencies ipit) are discussed). 

(ii) Periodic time dependence Lp{t) — cos{ujt). In that case, diffusion in the energy space is suppressed by quantum 
interference, leading to the phenomenon of dynamic localizatiorJ^^L^. We do not discuss the periodic problem in 
the present paper. 

We further specify to the case of linear motion along the trajectory (p{t) = vt and replace the parameter C(0) by a 
more convenient for our present discussion dimensionless parameter a defined as 



«=T^^(OK-TTf^') ■ (3) 
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Depending on whether a is much smaller or much larger than one, the transitions between levels may be described 
either as Landau-Zener transitions between the neighboring levels or as transitions in the continuum spectrum ac- 
cording to the linear-response Kubo formula. In both limits, the quantum-mechanical state effectively experiences a 
diffusion in energy, so that the energy drift over a large time T is given by 



[£;(T) -£;(0)]2 ^ 2i:>rA^ (4) 

The dimensionless diffusion coefficient D depends on a. A remarkable result of Wilkinson^ is that in the case of 
unitary random-matrix ensemble, in both limits of large and small a, the energy-diffusion coefficient is given by the 
same expression 

D{a) = a. (5) 

Recently, this problem has been analyzed further by one of the authorgSiiS with the use of a cr-model constructed 
from the Keldysh-integral averaging over random matrices H{t). The cr-model formulation of Refs. 0, ^3 contains, 
in principle, full information about the function D{a), and its diagrammatic expansion allows us to compute the 
perturbative series for D{a) in the limit of large a. For the driven Gaussian orthogonal matrices, the one-loop 

correction has the form^'^^: D{a) — a{l + S^\~^a~^^^ -I- . . . ), where — r(I/3) 3"^/^. For the unitary ensemble, 

the number of loops in the diagrams for calculating D{a) must be even, which corresponds to expanding in powers 
of 



where we took into account that, according to Sec. I VII the expansion parameter is TT~'^a~'^/^. In Refs. IglandT^ it 
has been found that ^2 = 0, and it has been further conjectured that all higher-order perturbative terms also vanish. 
In the present paper we shall verify that this conjecture holds up to the four-loop order: the value of d/^ obtained by 
numerical evaluation vanishes with very high accuracy {\d4\ < 3 x lO"**) which is a strong argument in favor of 

di = 0. (7) 

In the process of developing the diagrammatic expansion for D(a) we observe a remarkable cancellation of dia- 
grammatic vertices of order higher than four. This cancellation is proven to all orders in the perturbation theor y by 
combinatoric means. Therefore we find that, at the level of perturbative series, the non-linear cr-model of Refs. lOllOl 
is exactly equivalent to a matrix (j)'^-type theory. We show that the resulting theory can be obtained from the initial 
CT-model by applying a transformation analogous to the Dyson-Malee v^^i^'^ parameterization of quantum spin oper- 
ators. The equivalence of the cr-model to a kind of a (/)^-type theory may be a more important result than just a 
tool for calculating higher-order corrections in (jnj: in particular, this equivalence between two theories may have its 
counterparts for other types of non-linear cr-models, such as the one describing the diffusion of a quantum particle in 
a disordered mediaiiiiii^. We perform one straightforward generalization of this diagram cancellation to the case of 
arbitrary time-dependence of the control parameter ip{t). 

The rest of the paper is organized as follows. In Section ^ we describe the cr-model of Refs. Hinl In Section mH 
we develop the rules of the diagrammatic expansion. Further in Section IIVI we apply these rules for constructing 
the perturbative series for D{a). In Section we prove the cancellation of the vertices of order higher than four 
in the "rational" parameterization. We further formulate the new diagrammatic rules and the corresponding (j)^ 
theory. In the following Section IVII we apply the derived equivalence to explicitly write down the corrections to the 
diffusion coefficient . The numerical evaluation of the diagrams to the four-loop order is reported in Section IVIII 
In Section IVIIII we generalize our results to a more general situation of arbitrary dependence of ip{t). In Section llXI 
we demonstrate that our 0"* theory can be obtained from the initial a-model by transformation similar to the Dyson- 
Maleev transformation. We conclude our discussion in Section ^ Technically complicated details of the calculations 
are delegated to Appendices. 



II. KELDYSH SIGMA-MODEL 



The starting point of our analysis is the cr-model action derived in Ref. 9: for the unitary ensemble. The field 
variable is the operator Q which is the integral kernel in time domain with values in 2 x 2 matrices in retarded- 
advanced Keldysh space. Technically, we write Q as a 2 x 2 matrix Qtt> depending on time variables t and t' . The 
general expression for the action with arbitrary ip{t) has the form 

5[g] = ^Tr^Q-^!^Tr[^,Q]2. (8) 
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Here E is the energy operator with the matrix elements Eu' = iStt'df, the diagonal in time representation operator if 
is defined by the matrix elements tptt' — Stt>f(t'), and the trace is taken both over the Keldysh and time spaces. The 
commutator in the last term of Eq. © vanishes if (p{t) = const, and is generally nonzero for a time-dependent ip{t). 

To study the system's dynamics for the case of the linear perturbation ip{t) — vt, we prefer to simplify the notation 
by measuring time in the units of A^^. The resulting action in dimcnsionless units may be written as 

S[Q] = Se[Q] + S^M, (9a) 

where 

SE[Q]^7nTTK.t{EQ) = -^Ti-K f dt{di-d2)Qt,t, , (9b) 

2 J ti=t2=t 



5kin[Q] = -—TiK.t%Q] = — Ti-K J J dtdt' {t-t'YQu'Qt't- (9c) 

Here Trx.t and TrK denote the traces over the full Keldysh-time space and over the two-dimensional Keldysh space 
only, and a is the same dimensionless coupling constant as in (|3Jl and ©. 
The Q- matrix itself is subject to the constraint = 1, or, more precisely, 

j dt'Qu'Qvv = (^^f' ) • (10) 

Where we have introduced the "retarded" and "advanced" (5-functions ~ 5(ti —t2^e) with an infinitesimal shift 
£. This definition ensures the proper regularization^^ of the functional integral of exp(— 5'[(5]). 

The functional integration in Q is performed over an appropriate real submanifold of the complex manifold defined by 
the constraint H1(J|) (this procedure is standard in the sigma-model derivation both in Keldyshi^*ii and supersymmetrio^ 
formalisms). We take this manifold to be the orbit 

Q = U-^MJ (11) 

of the saddle-point solution 

- (f ) (12) 

under unitary rotations U . The integration measure [DQ] is then the usual invariant measure on the orbit. 

We want to stress that the present approach slightly differs from the scheme originally proposed in Refs. ISlOlHol 
The sigma-model derived in Ref. 0, having the same action ((S)), was formulated on a different manifold Qf = 
Up^U-^AUUp = Up^QUp, where Up is a non-Hermitian rotation which contains the knowledge about the fermionic 
distribution function [see Eq. IjAlOjl for an explicit form oiUp]- However, since calculating the energy-space diffusion 
coefficient D is essentially a single-particle problem, one can get rid of the distribution function F in the definition of 
the integration manifold and express the diffusion coefhcient D in terms of a certain correlation function of the fields 
Q. This refinement of the theory is presented in Appendix 1X1 

The quantity of our interest will be the diffuson X'^(i) defined by 

(Qui QVu) = J [DQ] e-'^'^^Qiti Qi^L = ^'5(ti -t2 + h- h)V,,.t, ih - U) , (13) 

where 

Q[fl=TTK{<7^Qt,t.) (14) 

are the off-diagonal elements of the Q-matrix in the Keldysh space. The form of the right-hand side in (|13|l follows 
from the invariance of the action S[Q] with respect to time translations (Qu' ^ Qt+st,t'+st) and with respect to the 
energy shift (Qu' Qtt'e"^*-*^* •*). Using the causality principles, one can further show (see Section UTTl and Appendix 
IXj) that the diffuson must have the form 

P^(t) =^?(t)exphP(77,i)], (15) 

where P{rj — 0,t)= 0. Further we may expand P{ri,t) in rj. The diffusion coefficient Dia) defined in |@J is given by 
the coefficient at P{ri,t) at r/^t (at small rj and at large t). The derivation is given in Appendix IXI Formally we may 
write 



1 

D(a) = — lim — -—rr 



V^it). (16) 

77=0 
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III. DIAGRAMMATIC EXPANSION 



To develop the diagrammatic technique, we exphcitly parameterize the unitary rotations U in Hll() by elements of 
the corresponding Lie algebra. This parameterization may be chosen in many different ways, and we write generally 



where the components of W are given by 



C = /"'(»'), (17) 



H'„. ^ I _l ] . (18) 



and f{W) may be represented as the series 

f{W) = 1 + W + W^/2 + c3W^ + C4W^ + ... (19) 

The unitarity of U implies bu' = b^^ and fiW)f{—W) = 1. We have chosen to write /^/^ in the definition 117|l so 
that the parameterization of Q contains the first power of /: 

= Kf{W) (20) 

(note that W anticommutes with A). 

Possible choices of the parameterization will be discussed in detail in Appendix IbI Using the parameterization H18|l . 
(|20|l . the integration measure [DQ] may be written as [Db Db]J f[b,b], where J/ is the Jacobian associated with the 
parameterization /. The explicit expression for the Jacobian in terms of the function / is given in Appendix IbI 

Then we do the algebra of substituting the parameterization (|12l) . H19() . 1201) into the action H9al) ~ (|5c) l to obtain 



^(>4)r;, l^ , e(>4) 

u [u, uj — u ■ - [u, f J -r »- 

where 



S[b, b] = 5(2) [6, b] + 4^') [6, b] + [5, 6], (21) 



^^'^ ^] = i // '^^i ^^^2 bi2 [{di + d2) + a{h - t2f] b2i (22) 



is the quadratic part of the action, and Sg"^^' [&, b] and S^hi^^ ^] higher-order terms (only even-order terms in 

6, b appear in the action). 

The propagator of the quadratic action S''-^) \b, b] is the bare diffuson: 

(6t+,/2,t-,/2^t'-.'/2,t'+.'/2)'°^ = l^iv - ^')2?r it - t') (23) 

with 

Vl°\t) = 9{t)exp[-a'n^t]. (24) 

Now we build up the diagrammatic expansion: by expanding the higher-order terms 4"^^ ^] ^^'^ ^ktn'^ ^] 
generate vertices which we then connect by the propagators (|23() of the quadratic theory. Each propagator (|23() will 
be graphically represented as a band with two incoming legs (corresponding to the 6-end of the diffuson) and two 
outgoing legs (the 6-end of the diffuson), see Fig. ^ These two ends are not equivalent: the 6- vertex must contain 
later times than the b vertex (due to the retarded 9{t) factor in the diffuson 124|l '). Correspondingly, we will draw an 
arrow on the diffuson line pointing in the direction of decreasing time, from b to b. 

The nonlinear vertices arising from S^^"^ [b, b] have the form 

4'") [6, 5] = (-l)"+Vc2„Tr[(96)6(66)"-i], (25) 

where we use the shorthand notation (96) 12 = (9i + 92)612, and the trace is understood as the convolution in time 
variables. The vertices from S^^^ [b, b] are of the form 

S^n^b, b] = ajdh... dhn 612623 . . . 62„,i «?(*^ - ^j)'' (26) 
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(b) (c) 



FIG. 1: (a) Bare diffuson Vi°>{t). The arrow points in the direction of decreasing time, (b) Vertices arising from S^J^*\b,b]. 
Crossed diffuson denotes differentiation db in the vertex 1251 . (c) Vertices arising from S^^\b,b]. The vertex itseff is charac- 
terized by a graph connecting different times ti at this vertex. The numbers on the links of the graph specify the coefScients 
a<;' in 1^. 



where a[j^ are numerical coefficients expressed through the coefficients Ck of the expansion H19|l : 

4"^ = |^(-l)'--''+"c,_.C._,+2n. (27) 

Graphically, we order the diffusons at each vertex clockwise according to their appearance in the traces H25|l and H26|l. 

The average value of the physical observables is given by certain correlators of the original Q-field. When we 
compute them in terms of the b and b fields, those observables generate "external vertices". The vertices generated by 
S^f^'' [6, b] and S"^!^^^ [b, b] will be further called "internal vertices" . Finally, vertices generated by the Jacobian Jf [b, b] 
(unless it equals one) we shall call "Jacobian vertices" (see AppendixlElfor the explicit form of the Jacobian vertices). 

There are two obvious rules for constructing the diagrams. Firstly, the diffusons may be drawn on a planar figure 
without "twisting" (but intersections of different diffusons are allowed, see, e.g., Fig.|21;). Secondly, the diagram must 
not contain closed loops formed by internal lines. Such loops would immediately produce the factor 9{ti — <2)6'(^2 — 
ts) . . . 9{tn — ti) which always vanishes [note that 9{ti — ti) = due to the causality rulei*^-, and so loops of length 
one are not allowed either]. For example, the diagrams in Figs.|2tb identically vanish and should not be considered. 
An example of a non- vanishing two-loop diagram is shown in Fig. 

From the above constraints on the diagram we can prove that only even-loop diagrams contribute to the diffuson 
(|15|l . More generally, the following theorem holds: Let 6(2"+^) and 5'-^"+^^ denote the operators 



^(2n+l) ^ j i^^^^fj^^^^fj^^^^ ^ ^ ^ fj^^^^^ _ ^ ^^^^^ ^28) 

Then any nonzero average (fe*-"^-* . . . 6^"*=^ 6*^"^^ . . . b^"*-*) must contain equal number of b- and 6-operators and its 
diagrammatic expansion contains only diagrams with even number of loops. 

Furthermore, from the same principle it is easy to prove that = if ii < t^ in all orders of the 

diagrammatic expansion. This proves the causality of the full diffuson i|13|l : !?,,(< < 0) = 0. 

In the main body of the paper, we shall only use the rational parameterization: 

The Jacobian for the rational parameterization is knowniiS to be equal to 1. (In fact, such a parameterization is 
not unique. As shown in Appendix^ the class of parameterizations with unit Jacobian is given by a one-parameter 
family, with the rational parameterization being a particular representative.) 



FIG. 2: Examples of diagrams. Diagrams (a) and (b) trivially vanish. Panel (c) shows a non-vanishing two-loop diagram 
(this diagram is one of those considered in Refs. 



(a) t + ri/2^ — ^t' + r]/2 



(b) 



^3-0 — OO 



FIG. 3: (a) Irreducible self-energy block S,,(t — t'). (b) The diagrammatic expansion of the pro-diffuson {hb) denoted by the 
solid propagator, (c) The diagrammatic expansion of the full diffuson (Q^^^Q^~^)- The irreducible blocks differ from E due to 
the presence of external vertices (shown by empty circles) containing higher powers of b and b. 

In the rational parameterization H29|l . the coefficients of the series p9f) are given by Cfc — 2^^^ {k > 1) and the 
coefficients a^J''' in the nonlinear vertices H2B|I have the form 

4f = - (-l)*-J'+"2i-2»7r. (30) 
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IV. DIAGRAMMATIC SERIES FOR THE DIFFUSION COEFFICIENT 

In this section we apply the developed diagrammatic technique to calculating the full diffuson H13|l and to further 
extracting the diffusion coefficient (fTE|l . 
The observables in (|13|l have the form 

g(+) = Mi)-C3&(^)+C56(^)-..., 

Q(-) ^ 5(1)-C3&(3)+C56(^)-..., (31) 

where C3, C5, ... are the coefficients in the Taylor expansion of the parameterization H19|l . and the operators 6*^"^ and 
are defined in Eq. H28I) . Respectively, the average (Q'^-'Q'^-') contains averages of different powers of b and b. 
Let us first discuss the simplest of those averages: 

2 

(^t+7?/2,t-V2^t'-'j72,t'+V/2) = -^iV - V')'D*r,it - t')- (32) 
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Eq. H32() defines the "pro-difFuson" 'D*{t) which should not be confused with the full diffuson given by Eq. (|13|l . The 
pro-difFuson 'D*(t) may be written in the perturbative series as shown in Fig.|3|3. Defining the irreducible "self-energy" 
block Ej, (t) (shown pictorially in Fig. by the condition that it cannot be split in two pieces by cutting any one 
diffuson, we can resum the diagrammatic series for the pro-diffuson as 

[p;M]-' = [pf (33) 

Here we adopted the frequency representation, with 

/ + 00 1 
e^-'vl°\u;)dt= ^ . (34) 

The role of the self-energy S,,(aj) is shifting the pole of the pro-difFuson 'D*{uj). 

One can see that the bare diffuson (|34(l in the frequency representation looks exactly as the diffuson propagator 
l/(— -|- Dq^) for a particle in a disordered medium. However, the analogy between these problems does not extend 
beyond the formal coincidence of the bare propagators. Indeed, for a particle in a random potential, the frequency 
u) and momentum q are "decoupled" from each other: each diffuson in a diagram has the same frequency, while the 
choice of their momenta qi is dictated by the momentum conservation law in the vertices. Contrary, for the dynamic 
problem, the times rji and tj of the internal diffusons 'Dl^.\ti) which constitute the self-energy block S,j(i) are linearly 
related to each other (see Section IVll for their explicit form). Therefore, in the frequency representation, the internal 
diffusons in I],j(aj) will enter with different frequencies Ui coupled to the time arguments r]j. Thus, the frequency 
representation 1)34(1 is not a way to calculate S,,(t), but is a convenient tool for summing the geometric series (|33f) . 

Note that the full diffuson 2?,,(t) differs from the pro-diffuson 'D*{t) by higher-power averages (6'"' • (see 
Fig-Ot)- We observe that all irreducible parts [both S^(t) and the irreducible blocks in Fig.O;] remain exponentially 
decaying with time even at ry = (in the frequency representation, they are non-singular functions of r/ and uj at 
77 and u! 0)^8,. Therefore we may sum all the diagrams to the form 

T. ^ ^ Z{ri,uj) 

T^ni.^) = —— 2 FTT' (^^^ 

—tuj + arj^ — 2j?)(w) 

where Z{ri,uj) is a regular function at ^ and w ^ 0. As we shall see from the further explicit calculations, 
S^(tj) ^ as 77 — > 0. Furthermore, the full diffuson T>j^{uj) remains unrenormalized at = (see Appendix IXI for a 
proof). Therefore, Z{ri = 0,cj) = 1. 

Thus the correction to the diffusion coefficient is determined by the leading term in the expansion of S,,(ti;) in 
and w: 

D{a) = a — lim — — — ^ = a — lim — [ S„(t) dt, (36) 

r;^0 77^ jj^O ry^ Jq 

where the lower limit of t-integration follows from the causality of S^(i) which can be proved analogously to the 
causality of the diffuson I?,,(t). 

The diagrams representing the irreducible block S^(t) may be classified in the number of loops L. From simple 
power counting (as shown in Ref. ^ the L-loop diagrams give a contribution to D{a) proportional to a^~^^^, cf. 
Eq. Therefore, the diagrammatic expansion is a series in the number of loops, and (since the number of loops L 
must be even) the expansion small parameter is a~^^^, up to some unknown number. A more accurate analysis of 
Sec. IVII shows that the actual small parameter is n~'^a~^^^. 

In the following section we shall see that certain different diagrams with the same number of loops L may partially 
cancel each other. 



V. CANCELING VERTICES OF ORDER HIGHER THAN FOUR 



In this section we show that, in the rational parameterization H29|l . in diagrams containing only internal vertices 
(e.g., in S^(t)), all vertices of order higher than four are cancelled in all orders of the diagram series. 

The outline of the calculation is as follows. The vertices with derivatives originating from 5'^-'*^ (Fig-Qt) generate 

terms dtD^^t). We transform those terms according to the equation of motion for the diffuson: 

dtV'^^\t) = 5{t) - aifv'^^\t). (37) 
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The resulting expressions have the form of (ti — tjY vertices in Fig.^ and may then be combined with the vertices 

originaUy generated by s'^^^^ \b, b] . 

When performed in the rational parameterization H29|) , this procedure leads to the cancellation of all vertices of order 
higher than four in all orders of the perturbation series. The calculation is somewhat technical with the combinatoric 
counting of diagrams and coefficients. This calculation is reported in detail in Appendix lO 

Of course, at the end of the above simplification procedure, the fourth-order vertices get modified from their original 
form. The final form of the fourth-order vertex is 

- TTQ; f — — 

Sl[b, ^] = — ^ dt2 dt3 dU ih ~ t2){h - U) bi2b23h^bii. (38) 

The resulting theory in terms of &- fields (which we will further call "6-theory" ) has the action 

SMb,b] = S'^^\b,b] + Sl[b,bi (39) 

where S^^^[b,b] is given by (|22() . 

We need to make two comments on the above derivation. 

First, up to now the correspondence is established only between diagrams containing internal vertices. The cor- 
respondence between physical observables (i.e., external vertices) will be discussed in Section |IXI In the further 
discussion, we employ the 6-theory for calculating S^(f) which contains only internal vertices. 

Second, when constructing the diagrams for S,((i) in the 6-theory, one of the (i2i-i — t2i) factors in Sl[b, b] always 
coincides with the external time difference rj. The integrals over times U converge^^ and therefore 

S^=o(t) = 0. (40) 

Obviously S,,(t) is an even and regular function of rj. Thus the Taylor expansion of S,j(i) in small 77 starts with 77^. 
The leading term in this expansion determines the correction to the diffusion coefficient, according to H36|) . We shall 
perform an explicit calculation in the next section. 



VI. DIAGRAMMATIC EXPANSION FOR D{a) 

In this section we classify the diagrams for I]^(t) in the 6-theory (|39|l and write down explicit integral expressions 
for the correction to the diffusion coefficient H36() . 

Consider any diagram (in the 6-theory) with L loops contributing to S,, (t) . As mentioned previously, such a diagram 
must contain no closed threads, which greatly restricts the number of allowed diagrams (in particular, the number 
of loops L must be even). Since all the vertices in the 5-theory have valency four, the number of diffusons in the 
diagram is 2L — 1. We enumerate those diffusons in two different sequences which we shall call the "left-hand" and 
"right-hand" numberings. Start at the "entry point" (external diffuson leading into the diagram) and go along the 
diffusons in two different ways: using the right-hand rule (following the right edges of the diffusons, i.e., turning right 
at every vertex) and using the left-hand rule (following the left edges of the diffusons and turning left at every vertex). 
The relative re-numbering of the diffusons in the two sequences defines a permutation a of (2L — 1) elements. In 
other words, let us number the diffusons following the right-hand route. Then reading off the diffuson numbers along 
the left-hand route produces the sequence of numbers [cr(l), cr(2),. . . , a{2L — 1)]. We shall further label the diagrams 
by such sequences (see Table QJ. An example of the 4- loop diagram corresponding to the permutation a = [3764215] 
(numbered 2 in the Table is shown in Fig.^ 

Note that not every permutation of 2L — 1 elements defines a diagram. We do not discuss here the combinatoric 
problem of enumerating all diagrams in all orders. For our modest purpose of calculating the four-loop correction to 
the diffusion coefficient, the enumeration of the diagrams may be easily done "by hand" . Recall that for L = 2 there 
is only one allowed diagram (shown in Fig. [21) corresponding to the permutation [3,2,1]. For L = 4, the number of 
allowed diagrams (and the corresponding permutations) is 21, see TableQ] Note that one of those diagrams (No. 1 in 
Table P) is reducible: it consists of two independent two- loop blocks. 

Given a diagram characterized by a certain permutation ct, let rji and ti denote the parameters of the diffusons 
'^Vi^ i'ti) involved in the diagram, with the diffusons enumerated along the right-hand route. For calculating the 
diagram, it is convenient to take ii,...,t2L-i as independent integration variables. The parameters rji may be 
expressed via ti as 

r,,=^ + Y,tj^ ^ = 77 + ^ ^^;k,. (41) 

j<i a(])<cr(i) 3 
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No. 




comments 


1 


[3214765] 


reducible 


2 


[3764215] 






3 


[3752164] 






4 


[3621754] 






5 


[3654721] 






6 


[4276315] 


= No. 


3 (T) 


7 


[4731652] 






8 


[4317625] 


= No. 


4 (LR/T) 


9 


[4726531] 






10 


[5274163] 






11 


[5417362] 


= No. 


3 (LR) 


12 


[5372641] 






13 


[6327514] 


= No. 


7 (T) 


14 


[6514732] 


= No. 


2 (LR) 


15 


[6251743] 


= No. 


3 (LR+T) 


16 


[6437251] 


= No. 


9 (LR+T) 


17 


[7614325] 


= No. 


5 (LR/T) 


18 


[7532614] 


= No. 


9 (T) 


19 


[7426153] 


= No. 


12 (LR/T) 


20 


[7361542] 


= No. 


9 (LR) 


21 


[7254361] 







TABLE I: The list of all four-loop diagrams for the self-energy part E. Each diagram is characterized by a permutation a, see 
text. The diagram No. 1 is reducible: it can be written as T,2'D^^^T,2, where E2 is the two-loop contribution. The other diagrams 
are irreducible. The diagrams related by the "left-right" (LR) or time-reversal (T) symmetries give the same contribution. 
LR-I-T means that two diagrams are related by the combined action of the two symmetries. LR/T indicate that two symmetries 
applied to a permutation give the same result. Graphical representation of the diagram No. 2 is given in Fig.0] 




FIG. 4: An example of the 4-loop diagram for the self-energy E. This diagram corresponds to the permutation a — [3764215] 
hsted as No. 2 in the Table |l| 



Here the first (second) sum contains aU times tj encountered before ti along the right-hand (left-hand, respectively) 
route. The last equality is the definition of the coefficients 7^^'*- The matrix 7^^^ contains only entries 0, 1, and — 1, 

and is antisymmetric: j^^"^ = Moreover, the elements j^^"^ with i > j are non-negative (either or 1), and 

those with i < j are non-positive (either or —1). 

Finally, calculating the self-energy block (i) in the 6-theory (|39|l and employing Eq. H36|) , we may write down the 
contribution of any given irreducible diagram (corresponding to a given permutation a) to the diffusion coefficient 
D{a): 

6^^^D{a) = -^a'-^/\ (42) 

-77- i/ 
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No. 






St. dev. 


2 


2 


0.066945 


0.000015 


3 


4 


-0.123205 


0.000012 


4 


2 


0.082358 


0.000018 


5 


2 


0.008234 


0.000005 


7 


2 


0.006030 


0.000007 


9 


4 


-0.006567 


0.000005 


10 


1 


0.206687 


0.000013 


12 


2 


-0.001646 


0.000006 


21 


1 


-0.011374 


0.000003 



TABLE II: Results of numerical evaluation of the diagrams from Table|I] M'°^' denote the multiplicities of the diagrams. K' 
is given by Eq. 1)43^ . The last column is the standard deviation of K^'^\ 



where 



Here P^'^\Ti, . . . ,T2l-i) and S'3'^'(ri, . . . ,T2l-i) are homogeneous polynomials of degrees 2L — 2 and three, respec- 



ts 

lively. They are defined as 



(2L-1 \ 2L-1 
n ^0 ^ (^^^^ 
i=i / i=i 

2L-1 

si''\T,,...,T2L^,)=J2vfT^, (44b) 

i=l 

fl^ = T.^^T,, (44c) 



where the coefficients ^["^^ are constructed from the permutation a according to Eq. H41|l . 



VII. NUMERICAL CALCULATION OF THE 4-LOOP DIAGRAMS 

The list of allowed four-loop diagrams (permutations) is given in Table n The first dia gram in the list is not 
irreducible: it consists of two blocks of two loops each and should not be included in the irreducible part S,((t). Out 
of the remaining 20 diagrams, some are related by symmetries and produce the same contributions 1)43(1 . Namely, there 
are two possible symmetries: the "left-right" (LR) reflection of the diffusons (corresponding to the transformation 
a ^ cr"^) and the time-reversal (T) symmetry (corresponding to replacing <j{i) ^ 2L — a{2L — i)). The diagrams 
related by the symmetries are identified in the above table. It remains to calculate the nine different diagrams and 
add them with the corresponding multiplicities M^''\ see Table HTl 

The antisymmetric matrices 7'-°^^ defined for nine permutations a from Table Ull are listed below, with the elements 
+ (— ) standing for 1 (—1) for brevity: 
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Here, the matrix ^'^"^^ denotes the matrix 7^°') for the permutation No. m from the Table UTI 

The resuhs of Monte Carlo numeric evaluation of the seven-fold integrals K^"") given by Eq. (|48|l are summarized 
in Table ITTl Performing summation with the multiplicities M^'^^ we get finally the estimate for the coefficient in 
the Taylor expansion 



di = -(7 ±7) X 10" 



(46) 



The numerical uncertainty indicated in Eq. H4ti|) corresponds to one standard deviation. Thus, we cannot distinguish 
di from zero and can estimate the upper bound for its absolute value as |d4| < 3 x lO""*^. 



VIII. ARBITRARY DEPENDENCE OF ip{t) 

In this section we discuss to what extent the results obtained above can be generalized to an arbitrary time 
dependence of the control parameter ip{t). 

We start by summarizing the modifications of the theory introduced by an arbitrary ^p{t). For a generic </3(i), the 
diffuson becomes a function of three times [cf. Eq. H13() ]: 



(47) 



To study the action ||SJ) one can develop the standard perturbation theory described in Section ITTll In terms of the 
&- fields, the action will have the form H21|) with the quadratic part {ipi = f{ti)) 



5(2) [6, 6] 



2A 



dti dt2 hi 



(9i + 92) + r((^i -^2)' 



&21 



(48) 



and infinite number of nonlinear terms S-'^[h,b]. In Eq. H48I) we introduced F ~ 7rAC(0). The bare diffuson defined 
as the propagator of S'(2) [b, b] is given bjsiSiSiiS^ 



(t, t') - B{t - t') exp \ -T / [^(r + 77/2) - ^{r - ii/2)fdT \ . 



(49) 
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It is remarkable that with such a modification of the theory one can still prove the cancellation of all internal 
vertices of the order higher than 4 in the rational parameterization. Indeed, the proof presented in Section [V] was 
based on the equation of motion H37(l for the diffuson, which is now replaced by an analogous equation 

a^Pf - 6{t - t') - r[^{t + r,/2) - ^{t - ri/2)fV^\tX): (50) 

and the combinatorial counting of coefficients which is insensitive to time dependence of ^{t). 
The resulting ^-theory has the action 

5 = - — j dtidt2dt'idti{ipi ~ ip2){^3 - ^i)bi2h2zh4bii, (51) 

which is a generalization of Eq. ()39|l for an arbitrary dependence of the control parameter ip{t). Thus, the diagrams 
for the pro-diffuson {hh): 

2A 

(&t+V2,i-V2^t'-'?72,t'+'?72) = ^ v')'D*{t,t') (52) 

with an arbitrary ip(t) are exactly the same as the diagrams in the linear case ip{t) — vt considered in the previous 
Sections (but with the new diffusons if^ l. 

The energy absorption rate is expressed through the full diffuson 'Djj{t, t') defined in terms of the field Q by Eq. l|T7|l. 
It differs from the pro-diffuson T>*(t,t') in the 6-theory since Q is a nonlinear function of b and b. In studying the 
linear perturbation this difference was irrelevant for the calculation of the diffusion coefficient (see Section rVT)) . For 
a generic perturbation (p{t), the difference between the full diffuson V^it^t') and the pro-diffuson 'D*{t,t') becomes 
important. In particular, the two-loop analysis of the quantum interference correction under the action of a harmonic 
perturbation-'^'' has shown that the average {bbb ■ b) (the diagram (c) inRef.[l3) has a contribution comparable to that 
of the average {b ■ b), both being negative and growing with time oc t. 



IX. RELATION TO THE DYSON-MALEEV TRANSFORMATION 



In this section we elucidate the meaning of the effective 6-theory H39|) and discuss its relation to the Dyson-Maleev 
transformation widely used in dealing with quantum spin ferromagnets. 

The action (|39|l has been obtained in the previous sections from an analysis of mutual cancellations of higher-order 
vertices in the perturbation theory for the rational parameterization. However, it turns out that the same action may 
be obtained directly from the initial cr-model action (j^ if one adopts the following parameterization of the Q-matrix 
in terms of the fields b and b: 

Q^(l-bb/2 b-bbb/A\ (33^ 



-1 



bb/2j 



Indeed, this parameterization respects the nonlinear constraint — 1, and the trivial algebra of substituting H53() 
into the initial cr-model action (|5J) immediately leads to the effective action (|39f) of the 6-theory. With the explicit 
parameterization (|53|l , external vertices (matrix elements of the Q-matrix) also become finite polynomials in b and b. 
We have checked that the two approaches (direct use of the parameterization ()53f) and the vertex cancellation by the 
technique developed in Appendix produce the same effective external vertices. 

Thus we come to an important conclusion about the cr-model @: The rational parameterization H29I) (which 
contains an infinite series of higher-order interaction vertices) after mutual cancellation of higher-order vertices in the 
diagrammatic expansion is perturbatively equivalent to the parameterization (|53() . 

In the rational parameterization [more generally, in any parameterization of the form H18|l - (|20|) ]. the matrices b and 
b are Hermitian conjugate of each other: b = b'' . The same remains therefore true for the 6-theory H39|l obtained from 
the rational parameterization. On the other hand, the parameterization (|53|) with b — b'' violates the hermiticity of 
the Q matrix. Nevertheless, our derivation indicates that this violation is inessential at the perturbative level and can 
be taken into account by a proper deformation of the integration contour over the elements of the Q matrix. 

The parameterization (|53|l is closely related to the famous Dyson-Malee\ii2ii^ parameterization for quantum spins. 
In that representation, the spin-5 operators are expressed by the boson creation and annihilation operators d^ and d 
as 



5+ = (25 - ata)a, = a\ S''=S-d^d. 



(54) 
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The Dyson-Maleev transformation conserves the spin commutation relations but violates the property (S*")^ = S'^ 
rendering the spin Hamiltonian manifestly non-Hermitian. This is the expense one has to pay for making the spin 
Hamiltonian a finite-order polynomial in boson operators [as opposed to an infinite series in the Holstein-Primakoff 
parameterization^^; note that the Holstein-Primakoff parameterization is analogous to the square- root-even parame- 
terization (jB16|l ] . The Dyson-Maleev transformation has proven to be the most convenient tool for studying spin wave 
interaction^^iSi: it reproduces all the perturbative results obtained with the Holstein-Primakoff parameterization in 
a much faster and compact way. 

A classical analogue of the Dyson-Maleev transformation was recently used by Kolokolov^^ in studying two- 
dimensional classical ferromagnets. To the best of our knowledge, there exists just one article by Gruzberg, Read and 
Sachdev^^. where the Dyson-Maleev parameterization was applied for the perturbative treatment of a cr-model (in the 
replica form). 

In the Dyson-Maleev representation, the Hilbert space of free bosons should be truncated in order for the operator 
to have a bounded spectrum. In Ref. |26L this truncation corresponds to integrating over a bounded region in 
bosonic variables. We expect that a similar constraint on the matrices b and b may be required in order to achieve a 
nonperturbative equivalence between the initial cr-model in the Q-representation and the 5-theory. This question is 
of importance for studying non-perturbative effects in cr-models, but goes beyond the scope of the present paper. 

X. DISCUSSION 

This work has appeared as a result of our attempt to prove the conjecture formulated in Ref. Isl that the Kubo 
formula for the energy absorption rate of a linearly driven unitary random Hamiltonian gives an exact result in the 
whole range of driving velocities. In the notation of the present paper, this conjecture implies that the relation 
D{a) = a is exact. 

Being unable to verify the conjecture nonperturbatively, we calculate the four-loop correction to the Kubo formula 
D{a) = a in the limit of large velocities of the driving field, a 3> 1. In the process of our derivation, we have 
refined the cr-model approach of Ref. |3 our improved method does not involve the fermionic distribution function 
and expresses the energy diffusion coefficient H16|l in terms of the full diffuson H13|l of the field theory l^. 

We have further proven that the resulting cr-model is perturbatively equivalent to the specific matrix (/)'*-theory 
H39|l . The lack of higher nonlinearities makes it possible to classify all four-loop diagrams and write down analytic 
expressions for them without resorting to computer symbolic computations. The final evaluation of emerging 7- 
dimensional integrals H43|) cannot be done analytically, and we calculate them numerically. We find that the coefficient 
di in the expansion © is indistinguishable from zero within the precision of our calculation, with its absolute value 
bounded by |d4| < 3 x 10""*. 

We believe that this conclusion acts in favor of the conjecture D{a) = a. 

This result may appear less surprising if we recall that the static (time-independent) unitary random-matrix en- 
semble is also known to possess some peculiar properties. In particular, its spectral statistics can be mapped onto 
the problem of one-dimensional noninteracting fermions; for a certain class of integrals over the unitary group the 
saddle-point approximation is exact (the Duistermaat-Heckman theoremS^, see Ref. |23for a discussion). However, in 
the time-dependent problem considered in the present work, we are not able to perform mapping onto free fermions, 
and the Duistermaat-Heckman theorem does not help to evaluate the functional integral of the Keldysh cr-model. The 
cr-model is nontrivial, and its diffuson 2?,j(i) is a complicated function of rj and t at arbitrary value of the coupling a. 
The diffuson T>,^{t) is renormalized from the simple diffusive form 124|l . and our result indicates that only its long-time 
asymptotics [which determines the energy diffusion coefficient (|16|l ] is free from perturbative corrections up to the 
order a~'^^^. 

The verification of the original nonperturbative conjecture D(a) = a remains an open question. 

As a byproduct of our analysis, we have rederived the Q-matrix parameterization (|53|l which is analogous to 
the Dyson-Maleev parameterization for spin operators^'. The parameterization H53(l is not specific to the Keldysh 
formalism, and can be applied to a wide class of unitary cr-models (e.g., to that describing diffusion of a particle in 
random media, both in the supersymmetric or replicEk^i approaches), leading to a considerable simplification of the 
perturbative expansion. We could not extend this parameterization to the orthogonal and symplectic cr-models, since 
those involve additional linear constraints on the Q-matrix, apparently incompatible with our parameterization. 
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APPENDIX A: WARD IDENTITY AND ENERGY ABSORPTION RATE 



In this Appendix we show that in the many-particle formulation the energy absorption rate under the action of an 
arbitrary perturbation (f{t) is expressed through the generalized diffuson (|47|l . In particular, for the linear perturbation 
(p{t) = vt, the diffusion coefRcient in the energy space may be extracted from the decay rate of the diffuson H13(l . As 
a byproduct of our discussion, we also prove that, at ry = 0, the diffuson reduces to the step function: 

P,=o(0 = 0{t) . (Al) 



1. Two approaches to energy diffusion 



The diffusion coefficient in the energy space may be defined in two different ways. 

In the original sigma-model derivation®, the states of the time-dependent Hamiltonian H(t) were occupied by 
noninteracting fermions. In such a multi-particle formulation, the step-like structure of the fermionic distribution 
bmction f{E), 

lim f{E) = 1, hm f{E) = 0, (A2) 

E^ — oo E^oo 



generates a spectral flow of fermions from low to high energies leading to the increase of the total energy of the system 
with time. At large time and energy scales, this many-particle process may be described by a diffusion equation on 
the distribution function f{E). The corresponding diffusion coefficient expressed as DA^ (where A is the average 
interlevel spacing and D is dimensionless) translates into the energy pumping rate W = DA^ (the density of states 
is A-i). 

On the other hand, the same diffusion process may be observed in a single-particle problem. In the single-particle 
quantum mechanics Q), the energy of the system E{t) defined with the instantaneous Hamiltonian 

E{t) = {^{t)\H{m{t)) (A3) 

diffuses with time as described by (0}. One can expect that in the process of time evolution of H{t), the relative phases 
of the wave function components corresponding to widely separated energies become uncorrelated, and therefore the 
multi-particle diffusion evolution of f{E) and the single-particle diffusion process give equivalent definitions of the 
diffusion coefficient D. The above reasoning is justified for a non-periodic evolution of H{t) (for example, for the 
linear evolution of the parameter ip{t) = vt); for a periodic perturbation studied in Ref. 8, the phase correlations 
become important which leads to dynamic localization. Nevertheless, we show below that the energy pumping rate 
W in the multi-particle problem may be expressed in terms of the single-particle diffuson for a rather general time 
dependence (p{t). 



2. Ward identity 

A helpful tool for our further derivation are the Ward identities generated by rotations of the integration variables 
Q in the functional integral of the sigma-model. 
Consider the functional 

Il[V] = J V-^QVe-^^^'''^^^[DQ], (A4) 

where V is an arbitrary matrix in the time and Keldysh spaces (not necessarily unitary), and the action is given by 
Eq. © . The integration is performed over the matrices Q of the form 1)11(1 , 1)12(1 . If the rotation by the matrix V is 
local (with Vtt' — > 5tt' at ]t\, \t'\ — > oo), it may be compensated by changing the integration variable Q ^ V~^QV in 
Eq. l|A4p producing no anomalous contribution, and therefore H[V] is independent of such rotations V. 

The multi-particle approach described in the previous subsection may also be introduced with the same rotated 
path integral (|A4p . but with the anomalous rotation matrix V (involving the distribution function f{E))^-. This 
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anomalous rotation will be discussed in detail below in the subsequent subsection, and in this subsection we derive 
the Ward identity generated by infinitesimal non-anomalous rotations V. 

Expanding V = 1 + A + . . . and taking the variation of Il[V] with respect to the infinitesimal generator A we get 

SU^[{Q),A] + ^{QTri[Q,E]A))+n^C{O){QTr{[ip,Q0Q]A))^O. (A5) 

In deriving the second term with [Q,E] we performed a cyclic permutation under the trace, which is equivalent to 
integration by parts and omitting the resulting boundary terms. This procedure is justified for local Au'- 

Various components of Eq. (jA5|) associated with different elements of the matrix A give the set of Ward identities 
related to the t^-invariance of n[t/]. Of particular importance is its off-diagonal (Keldysh) component n(+) generated 
by the Keldysh element 

St.tJuj, + St,-t,^u^t,{d3 + di)Vt,^t,(^^^±^,*-^^^ J d^^^{Q^^{Q^,Q,^f-^) = Q, (A6) 

where the diffuson with three times Vnit, t') is defined in Eq. I|47|) for an arbitrary time dependence ip(t), and we have 
used that, due to causality, {Qtt') = A = a^Su' ■ 



3. Causality of the diffuson 

The diffuson 'Djj{t,t') may be calculated by summing the diagrammatic series as described in Section IlIII Every 
line in the diagram is the bare retarded diffuson H49() . and therefore the full diffuson i') also equals zero when 
t < t'. 

Furthermore, using the Ward identity (|A6I) at ?7 = 0, we immediately find that the diffuson t') reduces to the 
step function oft — t' . Indeed, setting — nullifies the last term yielding dt'T>jj=o{t,t') = —6{t — t'). Integrating 
this equation and using the causality of the diffuson {'D,j{t,t') = for t < t'), we obtain I?,,=o(i,i') = (^{t — t'). 

For our discussion in the next subsection, in order to regularize the diffusion process at < ^ — oo, we consider the 
situation where the evolution of the Hamiltonian H{t) switches on at a certain time moment to- This can be modelled 
by a time dependence of ip{t) such that it remains constant at earlier times ip(t < to) — const. For constant (p{t), the 
last term in ljA6p vanishes and (similarly to the case 77 = 0) we obtain 

^V,it,t') + 6it-t')^0 for t'<to-^. (A7) 
We can integrate this equation in two domains: 

V^{t, t') = e{t - t') for t, t' <ta- -y, (A8) 

and 

hi 

V^{t, t') ^ Vr,{t, -00) for t' <to- Y 

The first property (jASjl guarantees that the diffuson "switches on" only after the moment to when the evolution of 
the H{t) starts. The second property (|A9() states that the diffuson with the initial time before tg does not actually 
depends on this initial time, which allows us to define the diffuson originating at i = —00. 



4. Distribution function and energy absorption rate 



To prove the relation between the multi-particle and single-particle definitions of the energy diffusion (see the 
first subsection of this Appendix), we consider the multi-particle kinetic problem in which the system is initially 
in a stationary state characterized by an arbitrary fermionic distribution function f^^\E). The evolution of the 
Hamiltonian starts at a time moment which is described by (p{t < to) = const. Within the Keldysh formalismiSi^, 
occupation of states by non-interacting fermions may be taken into account by rotating the matrix Q by the upper- 
triangular matrix^ 

iVp,o,)u' = r"' ^"'1 , (AlO) 
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where F^^} is the Fourier transform (in t — t') of F(*^)(£') = 1 — 2f^'^\E). The evolution of the distribution function 
with time may be read off as 

F = in(+)[y^(o)]. (All) 

Note that Vp{o) is not a local rotation (it does not vanish at t ^ ^ co), and therefore F undergoes a non-trivial time 
evolution. 

We may calculate the time evolution of F in the same procedure as the derivation of the Ward identity ljA6|l . The 
only difference is that the second term in Eq. (jA5|) originating from S'e[1^~"'^Q^] vanishes for V — Vpio). This can be 
most easily seen in the energy representation, where F'^'^^ is diagonal and evidently commutes with E. Then, applying 
the Ward identity (|A6|I one arrives at the result 



nW[t/^(0,]i2 = -2j di3di4i^,f *3'^*i-*..U-t3(-93+a4)A,-*.(^4^,^4^) , (A12) 



which translates into 



Ft+^/2.t-^/2 = V^{t,-^)F<-°\7^). (A13) 



This equation proves that the diffuson X'^(i, t') formally defined as a {QQ) correlation function in the field theory © 
indeed plays the role of the evolution kernel for the distribution function F(t, E). 

The energy absorption rate can be expressed in terms of the distribution function asSiS 

7rz 

W{t) = -- lim at9„F,+^/2,t-„/2- (A14) 

Using Eq. (jA13|) and the asymptotics F^'^^ (77) ^ 1/(17777) at — > which follows from the fermionic boundary conditions 
(|X2)l . we get 



1 d 92 
2A at 9^ 



Vr^{t,-oo). (A15) 

i)=0 



If we consider the situation of a linear time evolution if ~ vt switched on at a certain time moment the 
linearly growing contribution to the three-time diffuson I?,,(t, —00) equals that of the two-time diffuson 2?,,(t — to) of 
the translationally invariant theory H13|l . This can be easily seen from the diagrammatic expansion of the diffuson 
I?^(t, —00): the characteristic decay time of diagrams is t* = A~^Q!~^/'^, and switching on at time to plays the role of 
a soft cut-off for the vertices (jHU- Therefore, at time scales t — to ^ t* , the details of the switching-on process become 
unimportant, and we may replace T>jj{t, —00) in (jA15(l by 'Dn{t — to) . This proves the equivalence of our single-particle 
definition H16I) of the diffusion coefficient to the earlier multi-particle approach of Refs. I8l9lll0l 



APPENDIX B: PARAMETERIZATIONS OF THE Q MATRIX 



1. The Jacobian 



In this Appendix we discuss different parameterizations of matrices Q by elements W of the corresponding tangent 
space, according to (|20|l . With such parameterizations, the integration over the non-linear manifold of matrices Q 
reduces to that over the linear space of W (i.e., over the fields bu' and bu')- 

J [DQ] = J e-5'^(^)[L»iy], (Bl) 

where J/[W] = exp[— S'j(W)] is the Jacobian depending on the function f{W) in (|20|l . 

The integration measure [DQ] may be defined from the invariant metric on the space of Q matrices: 

(dQ, dQ) = - Tr dQ dQ. (B2) 

Similarly, the integration measure [Z?M^] = [Db Db] may be defined with the metric 



{dW, dW) = -Ti dWdW ^2Tt db db. 



(B3) 
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The Jacobian exp[— S'j(M^)] is the square root of the determinant of the metric {dQ, dQ) relative to {dW, dW). Further 
in the Appendix, we compute explicitly the determinant action Sj{W) for any given parameterization f(W) and find 
those parameterizations which have the unit Jacobian {S,j{W) — 0). 

The function f{W) satisfies the nonlinear constraint: 

f{W)fi-W) = 1. (B4) 

This constraint is solved by the condition that In f{W) is an odd function of W. 
The function f{W) must be regular at = 0, with the Taylor expansion 

f{W) = l + ciW + C2W^ + . .. (B5) 

Without loss of generality, we can normalize W by the condition ci = 1. Then necessarily C2 — 1/2. Other coefficients 
have a certain freedom. For example, one can chose C2n+i arbitrarily, then C2n+2 are uniquely determined by the 
constraint (|B4|I . We naturally define cq = 1. 

The quadratic form {dQ, dQ) may now, using the cyclic trace property, be written as 

{dQ, dQ) = - Tr df{W) df{-W) = ^ a„,,„, Ti{W''' dW VK"^ dW) (B6) 

with symmetric coefficients a„j.„2 — ctn2,ni- Using the tensor-product notation, we can write this as 

{dQ, dQ) = {dW ® dW, R), (B7) 

where 

R= 5] ^W^"^ (B8) 

"1 ,"2 

and (•, •) in the right-hand side of ((B7)| is understood as {A^B,C^D) ^ Tr{ACBD). The Jacobian equals (det R)^/^ 
and can be generated as the "Jacobian action" Sj{W): 

Sj = -^TilnK. (B9) 

The final step of the derivation is to express R in terms of /. Since R contains only commuting matrices W, it is 
convenient to replace R by the function R{x, y) defined as 

R{x,y) = c^rn.n.x^'y"', (BIO) 

ni,n2 

SO that 

R = R{W (^l,l(g)W). (Bll) 

By inspection. 



rii n2 



a. 



"1 ,"2 



fci=0 k2=0 

In terms of generating fimctions f{x) and R{x,y), this can be written as 



R{x,y) = - ^(^)-^(^) /(-)-/(-^) 



X — y X — y 



iM _ I f{y) 
fiv) V /(^) 

x-y 



(B13) 



Now we can write down the explicit formula for Sj{W) in terms of f{W). For a given f{W), we may expand 



= E /?n..„2^"^y"% (B14) 

X — y ^-^ 

ni,n2 

and then 

Sj{W) = -J2 f3nun2Tr{W^')Tr{W"'). (B15) 
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2. Potentially useful parameterizations 

For doing perturbative calculation in the sigma-model one has to chose a certain parameterization. We want to 
mention four possibilities: 

• Square-root-even parameterization (in this parameterization, f{W) does not contain odd powers of W higher 
than one): 



f{W) = y/l + W^ + W. (B16) 

This parameterization is widely used in literature since perturbative expansion in this parameterization directly 
corresponds'^^ to the standard diagrammatic technique with cross averaging over disorder-23. The Jacobian in 
the square-root-even parameterization in not equal to 1. 

• Exponential parameterization: 

/(ly) = exp(VK). (B17) 

This form was used by Efetov to construct his famous parameterization^^ of the integration manifold for the 
zero-dimensional supersymmetric sigma-model, that opened a way to exact evaluation of the two-level correlation 
function. The Jacobian in the exponential parameterization is not equal to 1. 

• Rational parameterization: 

fm ^ (BIS, 

This parameterization was suggested by EfetOAii^ in studying the supersymmetric sigma-model. The rational 
parameterization is frequently used for its Jacobian is equal to l^iiS. 

• Square-root- odd parameterization (in this parameterization, f{W) does not contain even powers of W higher 
than two): 

/W = l + -^ + Wl + ^= I + ^ + ^ 1 ■ (B19) 

This parameterization was used in Refs. I29ll83l85l for the study of the level statistics of random Hamiltonians 
within the replica formalism. The square-root-odd parameterization is known to have Jacobian equal one. In 
addition to that, the choice of the parameterization IjBlQp renders the action of the zero-dimensional sigma-model 
S[Q] cx Tr AQ to be Gaussian in W. 

Below we shall prove that the class of parameterizations with unit Jacobian (e.g., Sj{W) = 0) consists of a one- 
parametric family (|B24|) . with the rational and square-root-odd parameterization being particular cases. 




3. Parameterizations with unit Jacobian 



Now we shall solve the problem of finding all parameterizations f{W) leading to the trivial Jacobian Sj{W) = 0. 
Some powers of W in the Jacobian action ljB15|l give zero traces: Tr(VF°) = due to the causality and Tr(M^") — 
for all odd n. Therefore we look for functions f{x) producing in ljB14|l only nonzero Pn-^^n^ with one of the indices 
zero or odd. 

Thus, a parameterization with unit Jacobian must satisfy the condition 



In^^^^^^^ ViW =27^^(2;) +Fi(y)+ ^2 (x,y) (B20) 

with antisymmetric F2{x, y) = F2{^x, ~y) = —F2{~x, y) = —F2{x, —y). By setting x = or y = we obtain 
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Finally, by symmetrizing the left-hand side of Eq. (jB20|) . we get rid of F2{x, y) and obtain a closed equation on f{x): 
[}{x) + f{-x)] - [f{y) + f{-y)] [f{x) + f{-x) - 2][,f{y) + f{^y) - 2] 



x'^ — y"^ x'^y'^ 



(B22) 



The latter equation may be simply solved as 



fix) + fi-x) = 2 + — (B23) 

with an arbitrary constant A. For f{x) this gives 



1-Af 



1 + 

/(^) = ^ -p —■ (B24) 

1 - A^ 

Thus the parameterization with the unit Jacobian form a one-parameter family. Out of the four examples of pa- 
rameterizations mentioned above, two belong to this family: at A = 1, the expression (|B24|) gives the rational 
parameterization, and at A = it gives the square-root- odd parameterization. The square-root- even and exponential 
parameterizations do not belong to this family and have non-trivial Jacobian contributions Sj{W). 

The main calculation of the paper (cancellation of higher-order vertices) is done in the rational parameterization. 
We could not obtain similar results in other parameterizations. 



APPENDIX C: CANCELATION OF HIGHER-ORDER VERTICES IN RATIONAL 

PARAMETERIZATION 



1. Combinatoric coefficients and diagrammatic rules 

The goal of this section is to show that all internal vertices of order higher than four cancel in the rational 
parameterization . 

In this Appendix, to simplify the figures, we shall pictorially denote the diffusons by single lines (with arrows), 
instead of double lines as in Figs. ^ and |21 The order of the arrows at any vertex remains important: it represents 
the order of 6-operators in the corresponding product Tr(666. . .b). In particular, arrows going to/from any vertex 
alternate between incoming and outgoing directions. Only vertices of even order are allowed. The diagram from 
Fig. 13; is again shown in Fig.[S^ in the new notation. The vertices may contain differentiations represented by a cross 
(Fig. ISJj) and [cpi — (pj)'^ terms (Fig. ^). When we transform a diffuson according to the equation of motion H37|l . 
there appear terms with the diffusons replaced by the 5-function. Such a replacement will be denoted further as the 
double-crossed diffuson, Fig. 01. 

Before turning to canceling vertices, we calculate the numerical prefactors at each diagram. For simplicity, we 
perform now the calculation only for diagrams representing corrections to the pro-diffuson H32(l defined as the (bb) 
propagator, with one incoming and one outgoing diffusons. Let No denote the number of diffusons, Ny the number 
of vertices, and 2ni be the vertex valencies. Then Njj = 1 + ^ and the number of loops in the diagram is 
L ^ Nd - Nv - I. 

Every diffuson in the diagram brings in an additional factor of (2/7r). In the rational parameterization 129|l . 

collecting the coefficients of the vertices S^f^^ [Eq. (|25|l ] and ■S'^^^'' [Eqs. H26() and H30() ] (and taking into account 
the n^-fold symmetry of vertices ()26f) ). the total numerical prefactor in front of the diagram for the pro-diffuson 

V* = (7r/2)(66) [Eq. ^] becomes {-l)^'^~^{2Tr)-^. With this prefactor, each vertex coming from S'g"^'' enters with 

the coefficient 1, and each vertex originating from 5'^j7j^'' enters with the coefficient aaij, where 

a., - (-1)-^+^ (CI) 

Finally, by rescaling the time in the units of a~^/^, the diagram acquires the overall factor a^^/^ (a^^^ for every 
vertex and from — 1 independent time integrations). Now we may factor out the common coefficient 

for all the diagrams of a given order L, and the remaining coefficients equal plus or minus one. Note that within a 
given order L, the relative sign alternates between diagrams with different number of vertices. 

The diagrammatic rules now may be formulated as follows. 
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FIG. 5: (a) The same diagram as in Fig. |5J; in the "thin-line" notation, (b) Differentiation dtD^ {t). (c) Vertex with 
{ifi — ipjY terms. At the vertex, the solid/dashed lines denote aij — ±1 respectively, cf. Eq. <CH . (d) Contraction of the 
diffuson: T>^r^\t) is replaced by 5{t). (e) An original vertex of the sigma- model action (four-valent vertex is shown as an 
example). The differentiations d and aijiipi — ipjf' terms are added together, (f) Graphical presentation of the diffuson 
equation of motion dtV':°\t) = 5{t) - rf'D'^°\t) (with a set to 1). 




FIG. 6: (a) Tree-like graph of contracted diffusons. (b) Loop of contracted diffusons. 

1. First we draw all the diagrams of a given order L allowed by the no-loop rule (see Section HlHl . The sign of the 
diagram is (—1)^°^^. 

2. In every vertex, we put the sum of all differentiations on outgoing diffusons (S's-generated) and the full graph 
of aij{(pi — ipj)"^ (S'kin-generated), see Fig.|S^. 

3. We expand every differentiation according to the diffuson equation of motion (Fig. |Sf). 

4. During this differentiation, some of the diffusons are replaced by S functions, which leads to merging or splitting 
some vertices (discussed in details below). Since every vertex has at most one differentiation on outgoing 
diffusons, the connected graph of "contracted" diffusons contain no more than one loop. We further distinguish 
two possibilities: 

(i) the connected graph of "contracted" diffusons is a tree (contains no loops). Fig. In this case, all the 
vertices of the graph are contracted to a single vertex ("root" of the tree; point A in Fig.|B^). As the result 
of summing such diagrams, the coefficients a.y at the "root" vertex get modified. 




FIG. 7: An example of a contraction tree and the corresponding polygon (room-painting) diagram. 




FIG. 8: (a) All admissible room partitions of the 8-wall room, (b) Every partition contributes the full graph of the "dis- 
tinguished" room A minus the external painted walls (and the whole graph is multiplied by ( — 1) to the number of internal 
walls) . 

(ii) the connected graph of "contracted" difFusons contains one loop, Fig.|Bl3. Then upon contraction, we obtain 
two disconnected vertices without any {^pi ~ ipj)"^ factors (vertices of this type would also be generated by 
the Jacobian action as described in Appendix IbI if the Jacobian were not equal to one). 

In the subsequent subsections we sum separately those two series of diagrams and show that in the case (i) 
vertices of order higher than four cancel, and in the case (ii) all such loop-contractions cancel each other. 
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A 



P(2)=l 




P(3)=-1 




+1 



P(4)=l 



FIG. 9: An example of a direct enumeration of p(2), p(3), p(4). 



Cancellations in tree-like contractions 



Consider now the case when the connected graph of contracted diffusons is a tree (Fig. This tree has one 

distinguished vertex (denote it A) which does not contain differentiation. All the {fi — LpjY terms of the resulting 
"contracted" vertex come from the vertex A. Now we represent the contracted vertex as a 2n-polygon, and the possible 
graph configurations leading to this vertex as dividing this polygons by several diagonals into smaller polygons (Fig.[71l. 
Each such diagonal (a "separator") should be thought of as a "wall" which represents a diffuson to be contracted. 
We "paint" such a wall from the side of the 6-operator (the side of differentiation). The "outside" walls of the 
big 2n-polygon are also originally painted in the alternating order (Fig. [71) . With this representation, the different 
contractible graphs generating a given vertex are in one-to-one correspondence with "admissible" partitions of 2n- 
polygons into smaller "rooms". A partition is called "admissible" if after painting all the walls in alternating order, 
no room has more than one painted separator. Obviously, with this rule, exactly one room of the partition has no 
painted separator. This room is called the "distinguished" room of the partition. As an illustration, in Fig. we 
show all the admissible partitions of the 8-wall room. From our discussion of the diagram rules and prefactors, each 
partition contributes to the contracted vertex the "full graph of the distinguished room" minus "painted walls of the 
distinguished room" with the relative sign of the parity number of rooms in the partition. Here the "full graph" 
means the full graph with the sign ^{—^y~^^'^{}Pi — ^ and the painted walls are subtracted when expanding the 
corresponding differentiations (see fig.|Hb for an example). 

Before we compute the total vertex as a sum of all room partitions, we define a combinatoric quantity p{n) . Consider 
a 2n-wall room with external walls painted in the alternating order, and one of the painted walls is "special" (say, 
it has a door). Let p{n) denote the "algebraic" number of all admissible partitions of the 2n-room such that the 
distinguished room contains special wall. By the "algebraic" number we mean that every partition comes with plus 
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FIG. 10: (a) A schematic picture of external blocks adjacent to the "distinguished" room A. (b) Combinatoric diagrams 
representing the external block configurations, (c) All possible external block configurations in computing p(4). 



or minus sign depending on the parity of the rooms in the partition (and the trivial partition with only one room 
comes with plus sign). For example, p{2) = 1, p{3) = — 1, p(4) = 1 (Fig. El, etc. 

One can prove by induction that p{n) — (—1)". Indeed, all partitions may be classified by "external" blocks 
(shaded in fig.llUbl adjacent to the distinguished room. Each of those external blocks contributes —p{ni) possibilities 
to subdivide it, where rii is the number of unpainted external walls of the block. All possible configurations of external 
blocks may be enumerated by dividing the total of n unpainted external walls into non-intersecting sequences of length 
rii. Thus, 

where the sum is taken over all possibilities to choose a set of non-intcrsccting blocks of adjacent elements from a 
sequence of total n elements (Fig. llOb ) . The only element excluded from the sum is the set where the whole sequence 
belongs to a single block. In Fig. llOb we illustrate the possible choices of such block sets for n = 4. Now we perform 
the induction step assuming that for all < n, p{ni) — (—1)"'. Then the same sum (|C2p may be rewritten as 
p{n) = ^(— l)"""*-!, where riadj is the number of adjacent pairs belonging to the same block. Finally, we note that we 
may parameterize block sets by independently choosing two neighboring elements either to belong or not to belong 
to the same block. This easily gives us the total sum. If all the sets were allowed, the sum would be zero. The only 
set which is not allowed is the one with all neighboring elements (total n — 1 pairs) belonging to the same block. This 
gives p{n) = (—1)", which completes the proof. 

Now we do the final step of the calculation: we calculate the renormalized coefficient dij with which the combination 
[ifii — ipj)'^ enters the contracted vertex. First, such terms with i adjacent to j through a painted external wall never 
enter (they are always subtracted by expanding differentiations). So we may assume that i and j belong to different 
painted external walls (let us call those walls Wi and Wj for future reference). 

A room partition gives a contribution to dij if and only if Wi and Wj both belong to the distinguished room. 
Similarly to our discussion above, for each room partition, we may consider external "blocks" complementary to the 
distinguished room. Now the walls Wi and Wj divide the unpainted walls into two sequences (define their lengths as 
m^^^ and m^^^), and the contributions from these two sequences factorize: 

a., = i-iy-^-' J2 Ei-pi^'^^)) E n(-^'("^i'^)) (C3) 

(1) (2) 
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FIG. 11: The painted walls Wi and wj divide the boundary into two pieces of lengths m*-^^ and m 




FIG. 12: The resulting four-valent vertex. 

(see Fig. Ill|l . The only difference from the previous calculation is that now the whole sequence m'^^^ (or m^^^) is 
allowed to be placed in one block. Therefore the sum gives zero unless both rrS^^ and m^^^ equal one. But this is 
possible only if the total number of unpainted external walls is n = ni^^^ + to'^^ = 2, i.e. for the four-valent vertex. 
This completes the proof that all higher-order vertices cancel in the process of tree-like contractions. 
The only surviving four-leg vertex is shown in Fig. 1121 with 

'Y^a^j{'^i - ipjf = {ip2 - ^Pzf + (<y3l - ipif - (<Pl - if^f - (<y52 - '-P'if = 2(^51 - <^2)(¥'3 " ¥'4)- (C4) 

The diagrammatic series with such a vertex may be considered as a perturbative expansion of the field theory H38|l . 
(123. 

3. Cancellation of loop-like contractions 

In this section we prove cancellation of vertices obtained in the process of contracting a connected graph of diffusons 
containing a closed loop (Fig. IHt*). After such a contraction, two vertices appear without any — (j)jY terms. 
This is the same type of vertices as those appearing from the Jacobian. We consider the rational parameterization 
which contains no Jacobian vertices, so the only vertices of this type are those produced by the closed-loop diffuson 
contractions. 

Here we make two important remarks. First, there is no internal structure of those vertices, and the only parameter 
of the vertex is its numerical prefactor. Second, it is sufficient to prove cancellation of the "loop part" of the diagram 
(Fig. I13|) containing only the closed loop of contracted diffusons, without any "branches" which do not depend on the 
internal structure of the loop part. 

So we consider loops of length ^ > 1 of contracted diffusons. If we go along the loop, the legs on the left and on 
the right sides of our path are collected into two new "contracted" vertices. We fix the valencies of these vertices to 
be 2n and 2m and count the total combinatoric factor p{n, m) as the number of ways to obtain these vertices from 
loops of different lengths. Taking into account the sign alternating with the total vertex number, we obtain 



oo 

p{n,m) = ^(-l)'p;(n,m), 
1=1 



(C5) 
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FIG. 13: The loop contraction of diffusons. The inner and outer boundaries of the shaded region in the right-hand side are 
contracted into two separate vertices. 




P2(l>2)=4 



FIG. 14: Direct computation of p(l, 2) = -pi(l, 2) + p2(2, 2) - p-i{l, 2) = 0. 

where pi (n, m) is the positive integer number counting the number of ways to create the two vertices of valencies 2n 
and 2m from the loop of length I. Note that from the left-right hand rule both n and m must be positive. Also, since 
the two-leg vertices are not allowed in the original action, piin^m) — for I > n + m. An example of calculating 
p(l, 2) is shown in Fig. 1111 

The combinatoric problem of computing pi{n,'m) may be related to computing another combinatoric quantity 
qi{n,m) defined as the number of ways to distribute n identical black and m identical white balls among I different 
boxes so that no box remains empty: 

pi{n,m) = nmjqi{n,m) (C6) 

(here factors n and m come from different ways to circularly re-label external legs, and I /I factors reflects equivalence 
of diagrams obtained by the circular permutation of vertices). The example of calculating qi{l, 2) is shown in Fig. 1151 
Finally, we use the generating- function method to compute qi{n, m). Consider the series 



fi{x, y) = {x + y + x^+xy + y^ + ...y, 



(C7) 
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qi(l,2)=l 

^ 1^ ^ ^ 

0«# •0# ##0 13(1'2)=3 

I I I I I I 

FIG. 15: Example of g;(n,m) definition: direct computation of 5;(1,2). 

where the expression in the brackets contains all terms x"^?/"" except for Ux = n.y = 0. Then the coefficient at x"y" 
in the Taylor expansion of f{x,y) equals qi{n,m). Now a straightforward calculation gives: 

1 



^'(-'^)n (i-.)(i^.) -v ' ^'''^ 

and then 

oo 

^ ji-lYfiix, y) = ln(l -x)+ ln(l - y) (C9) 

1=1 

(note the factorization!). Now the coefficient p{n^ m) equals nm times the coefficient at x"?/™ in the Taylor expansion 
of the expression (|C9|) . The latter however equals unless one of the two numbers n or m is zero. This finishes the 
proof: all vertices generated from loop contractions cancel. 
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